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ATP-dependent chromatin remodeling enzymes (CRE) are bio-molecular motors in eukaryotic 
cells. These are driven by a chemical fuel, namely, adenosine triphosphate (ATP). CREs actively 
participate in many cellular processes that require accessibility of specific segments of DNA which 
are packaged as chromatin. The basic unit of chromatin is a nucleosome where 146 bp '^ 50 nm of a 
double stranded DNA (dsDNA) is wrapped around a spool formed by histone proteins. The helical 
path of histone-DNA contact on a nucleosome is also called "footprint" . We investigate the mecha- 
nism of footprint traversal by a CRE that translocates along the dsDNA. Our two-state model of a 
CRE captures effectively two distinct chemical (or conformational) states in the mechano-chemical 
cycle of each ATP-dependent CRE. We calculate the mean time of traversal. Our predictions on 
the ATP-dependence of the mean traversal time can be tested by carrying out m-vitro experiments 
on mono-nucleosomes. 

PACS numbers: 87.16.ad,87.16.Sr,87.16.dj 
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I. INTRODUCTION 

A deoxyribonucleic acid (DNA) molecule is a linear 
heteropolymer whose monomeric sub units, called nu- 
cleotides, are denoted by the four letters A, T, C and 
G. The sequence of these nucleotides in a DNA molecule 
chemically encodes genetic informations. In the nucleus 
of an eukaryotic cell, DNA is stored in a hierarchically 
organized structure called chromatin [l|-[j] . The primary 
repeating unit of chromatin at the lowest level of the hi- 
erarchical structure is a nucleosome |5| . The cylindrically 
shaped core of each nucleosome consists of an octamer of 
histone proteins around which 146 base pairs (i.e.,^ 50 
nm) of the the double stranded DNA is wrapped about 
two turns (more precisely, 1.7 helical turns); the arrange- 
ment is reminiscent of wrapping of a thread around a 
spool. There are 14 equispaced sites, at intervals of 10 
base pairs (bp), on the surface of the cylindrical spool. 
Electrostatic attraction between these binding sites on 
the histone spool and the oppositely charged DNA seems 
to dominate the histone-DNA interactions which stabi- 
lize the nuclcosomes. Throughout this paper, the helical 
curve formed by the histone-DNA overlap will be called 
the "footprint". 

The DNA stores the genetic blueprint of an organism. 
If nucleosomes were static, segments of DNA buried in 
nucleosomes would not be accessible for various functions 
involving the corresponding gene |6|~[8|. However, in real- 
ity, nucleosomes are dynamic. Spontaneous dynamics of 
nucleosomes are usually consequences of thermal fluctua- 
tions whereas the active dynamic processes are driven by 
special purpose molecular machines called chromatin re- 
modeling enzymes (CRE) fuelled by ATP jol-fl^. Various 
aspects of chromatin dynamics has received some atten- 
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tion of theoretical modelers, including physicists, over the 
last few years [l5l - l27t . 

In general, "Chromatin remodelling" refers to a range 
of enzyme mediated structural transitions that occur dur- 
ing gene regulation in eukaryotic cells. To make DNA, 
which is wrapped around histone octamer, accessible for 
various DNA-dependent processes, it is always necessary 
to rearrange or mobilize the nucleosomes. In principle, 
there are at least four different ways in which a CRE can 
affect the nucleosomes [23: (i) sliding the histone oc- 
tamer, i.e., repositioning of the entire histone spool, on 
the dsDNA; (ii) exchange of one or more of the histone 
subunits of the spool with those in the surrounding solu- 
tion (also called replacement of histones) (iii) removal of 
one or more of the histone subunits of the spool, leaving 
the remaining subunits intact, and (iv) complete ejection 
of the whole histone octamer without replacement. Our 
theoretical work here is closely related to sliding. 

In the next section we describe a scenario in which 
either a CRE motor (or, other ATP-dependent motors 
that translocate along dsDNA) traverse the "footprint" . 
Because of the stochasticity of the underlying mechano- 
chemical kinetics, the footprint traversal time (FTT) is a 
fluctuating random variable. Extending an earlier model 
developed by Chou [23, we analytically calculate the 
mean FTT (MFTT) of the CRE motor. 

To our surprise, we found that the ATP-dependence of 
the various ATP-driven activities of CREs have not been 
studied systematically in the published literature. In par- 
ticular, we address the question: how does the MFTT of 
a CRE motor vary with the variation of the concentra- 
tion of ATP? This rate is not necessarily directly pro- 
portional to the rate of ATP hydrolysis by the CRE be- 
cause the mechanical sliding of the nucleosome need not 
be tightly coupled with the hydrolysis of ATP by the 
CRE. Therefore, we develope here an analytical theory 
predicting the ATP-dependence of the ATP-dependent 
footprint-traversal by CRE. We hope our result will stim- 



ulate systematic experimental investigations on the ATP- 
dependence of ATP-dependent CREs. 



II. CRE: PHENOMENOLOGY AND 
MOTIVATION FOR THIS WORK 

Chromatin is not a frozen static aggregate of DNA and 
proteins. Spontaneous thermal fluctuations can cause a 
transient unwrapping and rewrapping of the nucleosomal 
DNA from one end of the nucleosome spool; the cor- 
responding rates for an isolated single nucleosome are, 
typically, 4s~^ and 20 — 90s~^, respectively [SO]. In other 
words, once wrapped fully, the nucleosomal DNA remains 
in that state for about 250ms before unwrapping again 
spontaneously; however, it waits in the unwrapped state 
only for about 10 — 50ms before rewrapping again spon- 
taneously. Surprisingly, the accessibility of the nucleoso- 
mal DNA is only modestly affected if instead of a single 
nucleosome the experiment is repeated with an array of 
homogeneously distributed nucleosomes [Sll, |32| . More- 
over, folding of an array of nucleosomes makes the linker 
DNA about 50 times less accessible [Sll [S^]. Further- 
more, on the nucleosomal DNA, the farther is a site from 
the entry and exit points, the longer one has to wait to 
access it by the rarer spontaneous fluctuation of suffi- 
ciently large size [33|. Thus, nucleosomal DNA far from 
both the entry and exit sites is practically inaccessible by 
spontaneous thermal fluctuations. 

Can a nucleosome slide spontaneously by thermal fluc- 
tuations thereby exposing the nucleosomal DNA? In- 
terestingly, spontaneous repositioning of nucleosome on 
DNA strands is a well known phenomenon [SJI . How can 
one reconcile accessibility of nucleosomal DNA by such 
repositioning [34j with the difficulty of access by unwrap- 
ping from either end of the spool [33| ? If the DNA were to 
move unidirectionally along its own superhelical contour 
on the surface of the histone, at every step it would have 
to first transiently detach simultaneously from all the 14 
binding sites and then reattach at the same sites after 
its contour gets shifted by 10 bp (or multiples of 10 bp). 
But, the energy cost of the simultaneous detachment of 
the DNA from all the 14 binding sites is prohibitively 
large because the total energy of binding at the 14 sites 
is about IbksT [13, El [HI. 

But, why can't the cylindrical spool simply roll on the 
wrapped nucleosomal DNA thereby repositioning itself? 
If the nucleosome rolls by detaching DNA from one end 
of the spool, cannot it compensate this loss of binding 
energy by simulataneous attachment with a binding at 
the other end? If such an energy compensation were 
possible, detachment from only one binding site would be 
required at a time. But, the cylindrical spool has a finite 
size on which only a finite number (14) of binding sites 
for DNA are accomodated. Therefore, by rolling over 
the DNA the spool would not offer any vacant binding 
site to the DNA with which it can bind. This rolling 
mechanism would successfully lead to spontaneous sliding 



of the nucleosome only if the histone spool were infinite 
with an infinite sequence of binding sites for DNA on its 
surface 
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We now describe a plausible mechanism for sponta- 
neous sliding of a nucleosome (l6l - [l8j . In the process of 
normal "breathing", most often the spontaneously un- 
wrapped flap rewraps exactly to its original position on 
the histone surface. However, if the rewrapping of a un- 
wrapped flap takes place at a slightly displaced location 
on the histone spool a small bulge (or loop) of DNA forms 
on the surface of the histones. Since the successive bind- 
ing sites are separated by lObp, the length of the loop 
is quantized in the multiples of lObp [1^. Such a spon- 
taneously created DNA loop, can diffuse in an unbiased 
manner on the surface of the histone spool. In the begin- 
ning of each step DNA from one end of the loop detaches 
from the histone spool, but the consequent energy loss is 
made up by the attachment of DNA at the other end of 
the loop to the histone spool before the step is completed. 
Consequently, by this diffusive dynamics, the DNA loop 
can traverse the entire length of the 14 binding sites on 
the histone spool of a nucleosome which will manifest as 
sliding of the nucleosome by a length that is exactly equal 
to the length of DNA in the loop. The diffusing DNA 
bulge can be formed by a "twist" , rather than bending, 
of DNA |35l - [37| . Spontaneous sliding of a nucleosome, 
however, is too slow to support intranuclear processes 
which need access to nucleosomal DNA. 

It is now widely agreed that ATP-dependent active re- 
modeling of nucleosome can account for the fast sliding 
of nucleosomes. Nevertheless, bulging DNA loop is ex- 
pected to play a key role in the remodeling process [231 ■ 
Our model describes how a CRE motor can wedge itself 
at the fork between the histone spool and a transiently 
detached segment of dsDNA and, by exploiting the spon- 
taneously diffusing loop by an ATP-dependent ratchet- 
ing, traverse the footprint in a directed manner. Because 
of the intrinsic stochasticity of the mechano-chcmistry 
of the CRE and that of the diffusive motion of the DNA 
loop, the overall motion of the CRE is noisy and the time 
it takes to traverse the footprint is random. 

The main question we address in this paper is the fol- 
lowing: if < T > is the MFTT, what is the dependence of 
< T > on the concentration of ATP in the surrounding 
aqueous medium? To our knowledge, in the published 
literature neither systematic experiments nor any ana- 
lytical theory has addressed this question. In this paper, 
by extending Chou's model [2^ of CRE, we capture the 
role of ATP explicitly and derive an analytical expres- 
sion for the dependence of < T > on the concentration 
of ATP. 

A CRE may be regarded as a molecular motor where 
input energy is derived from ATP hydrolysis and the 
output is mechanical work. The directed movement of 
the CRE may be caused either by a power stroke or by 
a Brownian ratchet mechanism [2| . The kinetics of the 
ATP-dependent CRE motor is formulated in our model 
in terms of a set of master equations. The kinetic scheme 



can be interpreted in terms of both power stroke and 
Brownian ratchet mechanisms. 

One interesting question [2^ in the context of CRE is 
whether the CRE translocates along the DNA by mov- 
ing around the nucleosome, or whether the CRE anchors 
on the histone octamer and "pumps" DNA by pulling 
around the octamer. From the perspective of physicists, 
these two alternative scenarios can be viewed as merely a 
change of frame of reference- one is fixed with respect to 
the DNA whereas the other is fixed with respect to the 
CRE. Therefore, we describe the operation of the CRE 
with respect to a reference frame with respect to which 
the CRE translocates along the DNA; but, the model can 
be reformulated by a coordinate transformation so as to 
capture the alternative scenario where the CRE pumps 
the DNA. 



III. THE MODEL 
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FIG. 1: A schematic representation of the isolated nucleo- 
some. 

We model a mono nucleosome where a dsDNA is 
wrapped one-and-threc- fourth turn around a disc-shaped 
spool made of histone proteins (see FiglIJa)). Follow- 
ing Chou [2^, we consider the scenario where the CRE 
"wedges itself underneath the histone" . 

The sites of histone-DNA contact along the DNA chain 
is represented as a one-dimensional lattice. Therefore, 
the lattice constant is, typically, 10 bp ((see FiglTJb)). 
The total number n of lattice sites is equal to the total 
number of histone-DNA contact in a single nucleosome. 



A. Flap, loop and diffusive sliding of histone spool 

In this subsection we present a summary of Chou's 
ideas |29| which we need in the next subsection where 
we extend Chou's model. Here we consider the simple 
situation when no CRE is present and the kinetics of the 
system is governed solely by spontaneous thermal fluctu- 
ations. Because of these fluctuations, from either end of 
the histone-DNA contact region, small segments of DNA 
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FIG. 2: A schematic representation of the position of the ther- 
mally generated flap and it's diffusion along wrapped DNA. 



momentarily unwrap from the histone spool at a rate /c„. 
For energetic reasons, the most likely length of such a 
segment would be one lattice spacing, i.e., about lObp. 
Following Chou [23 , we call such unwrapped segments a 
"flap" . The rate of the reverse transition, in which re- 
binding of the DNA flap with the histone, takes place at 
a rate kb- 

A flap need not re-make the original histone-DNA con- 
tact. Instead; by pulling in an extra segment of the DNA, 
its next segment can bind with the last binding site on 
the histone spool, with rate a thereby forming what Chou 
[29t referred to as a "loop'. While located at either end 
of the lattice, a loop can revert to a flap at a rate /3. The 
rates a and /3 are well approximated by [29[ 
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where Ebend is the energy cost of bending the DNA into 
the shape of the loop. 

A loop can step forward or backward. In the absence 
of any CRE. the rates of the forward and backward stop- 
pings of the loop are equal (denoted by k) , provided the 
size of the loop Li^op remains unaltered (see FigU); in 
each forward step it unwraps one segment of DNA from 
the histone in the direction of its hop and re-wraps an- 
other equally long segment behind it. Therefore, one can 
approximate k by [29| 
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{kb + ku) 



(2) 



When a loop, after entering the lattice from one end, 
makes an eventual exit from the other end, it completes 
the "sliding" of the histone spool by a distance Lioop 
along the DNA in the opposite direction. Therefore, from 
the perspective of the sliding histone spool, its effective 
rate of hopping by a step of size Lioop along the dsDNA 
strand is the same as the rate p„ at which a DNA loop of 
length Lioop traverses the lattice of n sites from one end 
to the other. 

Suppose Vj (t) denotes the probability that the loop is 
located at j (0 < j < 1). Following Chou's arguments, 
based on master equations for Vjit), one gets [29| 
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In the absence of a CRE, the traversal of a DNA loop of 
length Lioop from left to right is as likely as that from 
right to left. Therefore, the histone spool can slide for- 
ward or backward, with equal rate p„, by a step of size 



Lioop- As we'll see in the next subsection, peeling off of 
the DNA from the histone spool by a CRE motor keeps 
decreasing the effective value of n which, in turn, in- 
creases the effective sliding rate p„ . 



B. Kinetics of CRE-driven directed sliding of 
histone spool 
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Next, we consider the effect of DNA loop diffusion on 
the ATP-dependcnt translocation kinetics of a CRE. The 
model and results presented in this subsection are exten- 
sions of Chou's work [29| by incorporating explicitly a 
Brownian ratchet mechanism for CRE motors. 

As in ref. [23 , we assume that the step size of the CRE 
motor is identical to the length of the thermally generated 
DNA loop. Therefore, the mechanical movements of the 
CRE motor can be described as that of a "particle" on 
the one-dimensional lattice on which the equi-spaced sites 
denote the histone-DNA contact points. We denote the 
position of the CRE motor on this lattice by the integers 
m. We now extend Chou's model [29| by exploiting a su- 
perficial similarity with the Garai-Chowdhury-Betterton 
(GCB) model (Sg for the Brownian ratchet mechanism 
of monomeric helicase motors. 

A DNA helicase unwinds a dsDNA and translocates 
along one of two strands. At any arbitrary instant of 
time, the configuration of the system looks very similar 
to that shown in Figlljb) except that the surface of the 
DNA spool and the dsDNA would be replaced by the 
two strands of the dsDNA itself. The lattice constant is 
Ibp in the case of a helicase whereas it is about 10 bp in 
Figlljb). In the Brownian ratchet mechanism, momen- 
tary local unwinding of a segment, typically, 1 bp long, 
takes place at the fork by spontaneous thermal fluctu- 
ation; the opportunistic advance of the helicase merely 
prevents closure of the segment. Similarly, in the Brown- 
ian ratchet mechanism of the CRE, the CRE is assumed 
to "wedge" itself just in front of the DNA-histone fork. 
The CRE motor can move forward only if the segment in 
front of it is unwrapped by thermal fluctuation. 



FIG. 4: A schematic representation of the position of the 
motor with two state of the model. 
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FIG. 3: A schematic representation of the transition of the 
CRE between its two states. 
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FIG. 5: A schematic representation of the Brownian ratchet 
mechanism. 



The mechano-chemical cycle of the CRE is captured 
in our model exactly the same way in which that of the 
helicase was formulated in the GCB model [3^. We as- 
sume that sequence of states in each mechano-chemical 
cycle of a CRE can be combined into two distinct groups 
which we label by the integers 1 and 2 (see FiglS]). The 
allowed transitions and the corresponding rate constants 
are shown in Fig|31 

The physical processes captured by these rate con- 
stants can be motivated by a comparison with the ab- 
stract Brownian ratchet mechanism, illustrated in FigIS] 
ATP hydrolysis by the CRE drives its transition from the 
state 1 to the state 2 at a rate a;2i . Let us assume that the 
motor experiences two different types of potentials in the 
states 1 and 2. Let us further assume that initially the pe- 
riodic potential, with asymmetric sawtooth-like period, is 
kept on for sometime and during this time the motor set- 
tles at a position that coincides with one of the minima of 
this potential. Now if this potential is switched off then 
the probability distribution of the position of the motor 
will spread as a syinmetric Gaussian. After sometime 
this Gaussian profile is broad enough to overlap with the 
next well (shadded region in the FiglS]), in addition to 
the original well. Now if the sawtooth potential is again 
switched on then, with a non-zero probability (that is 
proportional to the area of the shaded region) the motor 
will find itself in the next well. Our model accounts for 
this possibility with the transition associated with the 
rate constant w(2- There is also a finite probability that 
the particle stays back in its original well; this is captured 
by the transition with the rate constant uji2. 

The CRE motor would step forward at the rate w^fj if 



the next site in front is cleared. But if the next site is 
not cleared and it has to wait for the unwrapping of the 
DNA segment by thermal fluctuation. Consequently, its 
effective hoping rate 
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is reduced from the free hopping rate cj(2 by a factor that 
depends on both fc^ and kf,. 

When a diffusing loop reaches in front of the motor it 
momentarily creates a flap of two bond segments. Three 
different transitions are now possible (see Fig|6]): (i) the 
motor's position remains unaltered while the two open 
segments close, (ii) the motor moves forward by one step 
while one segment of the flap closes; (iii) the motor moves 
forward by two steps and the flap cannot close. The rate 
for the process (i) is ki,/2 irrespective of the "chemical" 
state of the motor. However, the rates of the processes 
(ii) and (iii) depend on whether the motor was in the 
"chemical" state 1 or 2. If the motor is in the state 2, the 
rate of the process (ii) is given by ((wfj)"^ + k^^)^^ and 
that of the process (iii) is given by ((01(2) "'^ + (1^21) "'^ + 
(^(2)"^)- Therefore, 
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with the normalization constant 
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where the symbols /o, /i and /2 are the probabilities of 
the processes (i), (ii) and (iii) above when the motor is 
in the "chemical" state 2. Similarly, 
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with the normalization constant 
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are the corresponding probabilities, when the motor is in 
the "chemical" state 1. 

Suppose, N is the maximum number of histone-DNA 
contacts possible in the nucleosome. Let m denote the 
instantaneous position of the motor, n is the distance be- 
tween the motor and the far end of histone-DNA contact. 
The master equations for the probabilities P{m, n, t) are 
as follows: 
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FIG. 6: Schematic representations of the possible transitions 
when the motor is in front of the open flap. In (a) and (b) 
the motor is in the chemical state 1 and 2, respectively. 



For n>N-f 1 



dPi(m,n)/dt ~ LL:i2P2{'m,n) — uj2iPi{m,n) 

+ w/2P2(™-l,ri-t-l) 

+ pN[Piim,n + l) + Pi{m,n-l) 

-2Pi(m,n)] (9) 



and 

dP2{m,n)/dt = uj2iPi{m,n) - uJi2P2{m,n) 
- oj(2P2{m,n) 
+ PN[P2im,n + l) + P2{m,n - I) 

-2P2{m,n)] (10) 

For n=N 

dPi{m,N)/dt = uji2P2{m,N)-uj2iPi{m,N) 
+ Uj{^P2{m~l,N + l) 
+ pN[Pi{m,N + 1) - Pi{m,N)] 
+ PN[hP2{m-l,N) 

+ giPi{m-l,N)] 
+ pN-igQPi{m,N-l) (11) 

and 

dP2{m,N)/dt = L02iPi{m,N)-uJi2P2{m,N) 
^ u{^P2{m,N) 

+ PN[P2{m,N + 1) - P2{m,N)] 
+ hPN-iP2{m,N -I) (12) 

For 3< n < N 
dPi{m,n)/dt = uJi2P2{'m,n) — L02iPi{m,n) 

+ uj(2P2irn-l,n + l) -pnPi{m,n) 
+ Pn[fiP2{m - l,n) + giPi{m - l,n)] 
+ Pn+i[f2P2{m-2,n + l) 
+ g2Piim-2,n+l)] 
+ gaPn-iPi{m,n-l) (13) 

and 

dP2{m,n)/dt = uj2iPi{m,n) - uJi2P2{m,n) 

- {ui{2+PN)P2{m,n) 

+ /oPn-iP2(m,n-l) (14) 

For n=2 
dPi{m,2)/dt = LOi2P2{m,2) - L02iPi{m,2) 

+ Co{2P2{m - l,i) - p2Pi{m,2) 
+ p2[giPi{m-l,2) + hP2{m-l,2)] 
+ P3[g2Piim-2,3) + f2P2im-2,3)] 

(15) 
and 

dP2{m,2)/dt = uj2iPi{m,2) - uJi2P2{m,2) 

- w{2-P2(m,2)-p2P2(m,2) 

(16) 
For n=l 
dPi{m,l)/dt = uJi2P2im,l)^uj2iPiim,l) 

+ Cd(2P2{m ~ 1,2) ~ kuPi{m,l) 
+ P2[l2P2{m - 2, 2) + g2Pi{m - 2, 2)] 

(17) 
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FIG. 7: A schematic representation of the position of the 
DNA-histone contact with two state of the model. 



and 

dP2{m, l)/dt = uj2iPiim, 1) - UJi2P2im, 1) - fc„P2(m, 1) 



(18) 



C. Footprint traversal time 

We define P^,„(i) = J2m ^m(™7 ": ^) ^^ the probability 
that the n histone-DNA contacts are intact at time t, 
irrespective of the position of the CRE motor. From 
equations (P|- (jl8p . summing over m, we get the following 
equations: For n>(N+l) 

dPi^n/dt = PN[Pl.n+l - 2Pi^n + Pl.n-l] + 

Wi2-f2,)i+l + '^12P2,n ~ '^2lPl,n (19) 

and 

dP2,n/dt = PN[P2,n+l - 2P2,n + P2.n-l] + 

^2lPl,n — '^12P2,7i ~ ^^12^2,71 (20) 

For n=N 

dPi^N/dt = -pnPi,n + 0J12P2.N + ^(2^2,^+1 

+90PN-lPl,N-l + giPNPl.N +PN 
Pi, N+1 + flPNP2.N - (^2lPl,N (21) 

and 

dP2,N/dt = -{lo(2+Pn)P2.N +0J2lPl,N -0Jl2P2,N 

+PNP2.N+1 + faPN-iP2,N-i (22) 

For 3< n < N 

dPl,n/dt = —PnPl,n + <^i2^2,n+l 

+U!i2P2,n — ^2lPl,n + 9QPn-lPl,n-l + 
giPnPl,n + g2Pn+lPl,n+l + flPnP2,n + 
hPn+lP2,n+l (23) 



and 



dP-Zji/dt = -{ui(2+PN)P2,n+^2lPl,n 
—^12^2,71 + /oPn-1^2,n-l 



For 11=1 

dPi.i/dt 
and 



-kuPi,i + '^(2^2,2 + 

f2P2P2,2 + g2P2Pl,2 + Wl2-P2,l 
^21^^1,1 



We define the survival probability S^^n{t) to be the 
probabihty that the CRE has not yet reached the far 
end of the footprint till time t, given that initially (at 
(24) t = 0) there were n intact contacts between the histone 
spool and the DNA on the footprint in front of the CRE 
motor. Obviously, Sf^^n{t) is the solution of the equations 
for Pf_i,„{t) with the initial condition ^^^^(O) = 1. 



(25) 
(26) 



dP2,l/dt = -kuP2,l + W2I A,l - Wl2-P2,l 

For n=2 

dPi,2/dt = -P2i^i,2 + ^{2^2,3 

+fflP2-Pl,2 + flP2P2,2 + g2P3.Pl,3 

+ /2P3-P2,3 + ^12^2,2 ^ ^21^1,2 (27) 

and 



Interestingly, the time-evolution of S^^n{t) can be re- 
cast as 



dSi^n/dt — a„ (S'l^n+i — Si^n) + a„ (5i^„_i — S'l.n) 

+rn{S2,n - 5i,„) (29) 



dP2,2/dt = -^12^2,2 + ^21 A,2 - P2J^2,2 
-'^12-^2,2 



dS2,n/dt — 6„ (>5'2,n+l — S2,n) + &„ (5'2,n-l " S2^n) 

+ Sn{Sl.n — 5*2, n) + Fn{Sl.n-l ~ 'S'2,npO) 



where the transition rates a„,&„,r„,s„ and Fn depend 



(28) on the value of n as follows: 
I 



For n > (iV + 1) 

Fn = '^{2' «rt = PN, a~ = PN, K = PN,K = PN, I'n = ^21, .S„ = UJ12 

For n=N 

Fn = Uj(2, a+ = goPN, a^ = PN,b+ = foPN, 6^ = PN, Tn = W21, S„ = (Wi2 + JiPn) 

For 3 < n < TV 

Fn = (^(2 + f2Pn),a,+ = goPn, a~ = g2Pn,bt = foPn, K = 0, r„ = a;21, S„ = {UJ12 -f flPn) 

For n = 2 

F2 = Cj{2 + f2P2,at = 0, 03^ = 52P2, bi^ = 0, 62" = 0, r2 = a;2i, S2 = ^12 + /1P2 

For n = 1 

Fi = 0,aj'" = 0,a^ = K^b^ =0,b^ ^ ku,ri = W2i,si = W12 



The master equations (|29)) - (p0)) together, effectively, cor- 
respond to the kinetic scheme shown in the Fig. [71 Using 
this scheme, the MFTT for the single CRE motor can 
be calculated analytically by extending the theoretical 
framework developed in rcf. [39| for calculating the mean 
first-passage time of random walks. 



by 



-f^jn — / ^fi,n\t)dt 



(31) 



Following Pury and Caceres [33], the MFTT is given Since S'^_„(cxd) = and 5^.„(0) = 1, integrating the equa- 



tions (|29l) and (pO|) with respect to t, we get 



Using (|43]) in (jM]) we, finally, get 



+ r„[T2,„-ri,„] (32) 

1 = ^Ti {T2,n+1 — T2^n) + b„ {T2^n~l ~ ^2^,1) 

+ SniTl,n - T2,n) + Fn(Ti^n-l " T2,n) (33) 



N 



Defining 



On = T2,n ^ r|_ „ 



equations (|5^ and (P5|) can be expressed as 

- 1 = a+Ai,„ - a;; Ai,„_i + r„J„ (35) 

-1 = 6+A2,„-&;;A2,„_i-s„(5„-F„(Ai^„^i+5„) (36) 
Now, in the special case 

Ai,„ = A2.„ = A„, (37) 

equations ([55| and ([M|) become 



1 = «« ^n - a., A„_i + r„(5„ 



(38) 



- 1 = 6+A„ - (fe- + K)A„_i - (s„ + K)<5„ (39) 

Next, multiplying equation ((38|) by (s„ + F„) and Eq. 
(EH) by r„, and then adding the resulting equations, we 
get 

~{rn + Sn + F.n) = {r„6+ + (.s„ + F„)a+ } A„ 

- {rn{b~ + Fn) + (s„ + F„)a-}A„_i 

(40) 

Eq. (POJ) can be re-written as 

- C„ = B„A„ - A„A„_i (41) 

where, 

C„ = (r„ + s„ + F„), 

-Bn = {^n&;t + (s„ + -F«)«7t}7 

^„ = {rn{bn+F,,) + {sn + F^)a-] (42) 

We can rewrite Eq. (|4ip as follows 

A„-i = fiA„ + ^ (43) 

For a fully wrapped histone, the MFTT td is given by 

N 

(44) 



td = ^^r^ 



td = Y^[Cn/An + Y^ C„+i/A„+, Jl B„+fe/A„+fe] (45) 



n=l 



fc=0 



Since it is not easy to get an intuitive feeling for the 
implications of the expression (|45p , we anaylyze its spe- 
cial simpler forms in some limiting cases. In the limit 
of extremely slow motor, i.e., ijJ\2 — >■ 0, as expected, the 
expression P5|) for the MFTT td diverges. 

For ensuring high-speed of the CRE motor, we need 

simultaneously w^j/^fc >> 1 ^-^^id W2i/ki, >> 1. If, for 

(34) simplicity, we make the additional assumption that 1/7(2 



is the slower of the two, i.e., W21 >> w^j, we have /o = 
/i = go = .91 — and /2 ~ 1, and 92 — 1- Hence, in this 
limit, Bn = for n < N and, therefore. 



TV 



td^J2 



n=l 



Cn 

A.r, 



N/ku 



(46) 



which is identical to the corresponding limiting value of 
td reported in ref. [23] . This is a consequence of the fact 
that in the limit of extremely fast motor, because of the 
assumption of very large value of 0121, the 2-state model 
reduces to an effectively 1-state model. We make a nu- 
merical estimate of td in this limit by computing an ap- 
proximate value of fcu . Defining 



K = kb/ku 



(47) 



as the flap binding constant, we can rewrite the equation 
dUl) as 



td 



NK 



(48) 



Range of typical values of K has been used earlier by 
Chou [2^. Using this range of values for K, one can 
estimate fc„, provided a typical value of kb is known. 
Therefore, we now estimate the typical numerical val- 
ues of kb following Schiessel and coworkers |16l4l8| . Sup- 
pose, Lo(— 50 nm) be the length of the DNA that wraps 
around the histone spool. Let L + dL he the contour 
length of the loop induced by spontaneous thermal fluc- 
tuations where (see Fig. [H]) i is the exposed arc length 
on the histone spool that was covered by the DNA seg- 
ment prior to the loop formation and dL is a small seg- 
ment of the linker dsDNA that has been pulled into the 
loop. 

We assume that the life time of a loop (r) is much 
shorter than the average time required to form a loop. 
Following Schiessel et al. pi5l - KL8| we write down the rate 
of loop formation as 



a 



Lgexpi-Ebend/jkBT)) 

tL' 
Comparing Eqs. ([1]) and P^ we obtain 

kb = —TT- 

tL 



(49) 



(50) 




FIG. 8: Top view of the histone octamer bound with DNA. 
Loop formation involving length dL of linker chain being in- 
corporated into the nucleosome, with leng th L' of the exposed 
surface, (adapted from Fig. 2 of ref . [lq|'). 
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FIG. 9: Normalized MFTT aN'^td plotted against the nor- 
malized motor speed (^12 /a for different values of K with 
UJ-12/0L = uj2i/ot = 0.5. 



Since t^^ — k characterizes the rate of unbiased (|ifFusion 
of the loop around the histone spool [3 ''' — "d where 
D is the corresponding diffusion constant. From Stokes- 
Einstein relation D ~ ksT/^, where C ~ 77L [1^ and 77 is 
the effective viscosity of the aqueous medium. Combining 
all the results and substituting these into Eq. (|50|) we 
finally obtain 



keT 
vLo{L'] 



(51) 



The estimation can be completed only if an estimate of 
L' is available. Following ref. [la] (Eq. (2a) of [la|), we 
get 



L c. (^)V6(dL/i?„ 



'"R, 



(52) 



where k is the bending elastic constant of the semi- 
flexiablc DNA chain, A is the adsorption energy per unit 
length and Rq is the radius of the histone spool. Us- 
ing the reasonable values quoted in ref. |16l4l8| . namely, 
i?o = 5nm,, K = 207.10piV — nm^,dL — 3AQnm and 
A = 5.92pN we obtain from Eq. §^ L' = 16A3nm. 
using this estimate of L , together with Ifc^T = 4.142 
pN-nm, Lq = 500A, and rj = 1 Centipoise, we obtain 
from Eq. ^ kb = 306877.4s~i. 

With the above estimated value of ki,, and A^ := 15 
from Eq. (gS]) we get the estimates td = 0.000024s for 
K = 0.5 and td = 0.0005s for K = 10. Such small 
values of td, estimated from EqHBJ arise from the fact 
that the approximate form (|46p is valid only in the limit 
of extremely fast motor. Therefore, this limiting formula 
provides only a lower bound and does not correspond to 
real CRE motors under physiological conditions. 

In FigH we plot the normalized MFTT tda/N as 
a function of the normalized motor speed (^(2/01 for 
Ldi2/a = uj2i/a = 0.5 and a few fixed values of the pa- 
rameter K. For any fixed value of K, the normalized 
MFTT decreases monotonically with the increase of the 
normalized motor speed and saturates to the value given 



by equation (^5)) in the limit 01(2/'^ ~^ ^^- Moreover, for 
a given value of ^(2/0^1 ^^ the flap binding constant K 
increases the MFTT increases. 

In FigfTUl td is plotted against oji2/a for (a) K = 10, 
uji2/a — 0.8, and (b) K = 10, a;i2/a — 0.1, each for a 
few distinct values of uj2i/ct- The MFTT decreases as 
u}2i/a increases. This is a consequence of the fact that 
W21 depends on the ATP concentration. For small iOi2/a 
reduces the amplitude of peeling time. 

In FiglTT]wc demonstrate that for large value of a;2i/a, 
which effectively speeds up the motor, reduces the mag- 
nitude of the MFTT. 

Although the qualitative trends of variations of td with 
w{n/a in our model is similar to that in Chou's model 
[29| , wide range of variation of td is possible in our model 
by controlling 0^21 which, in turn, can be controlled by 
the ATP concentration. 

In order to explore the dependence of td on the con- 
centration of ATP, we first assume that 



UJ21=UJ2MTP] 



(53) 



Assuming a typical value Wji = lO^M^^s^^, we have 
plotted the normalized MFTT against the ATP concen- 
tration for two different normalized values of the unhin- 
dreed motor speed keeping the other parameters fixed. 
With the increase of ATP concentration, the MFTT de- 
creases and, gradulally saturates. When ATP concentra- 
tion is sufficiently high, the step with rate constant LO21 
is no longer rate-limiting. We also find that, for a given 
ATP concentration, the higher is the value of 1^(2/'^ ^^e 
shorter is the MFTT td- 

The linear dependence of UJ21 on ATP concentration, 
as envisaged in (j53p . may be valid only at sufficiently 
low concentration of ATP. In general, W21 may follow 
the usual Michaelis-Menten equation for the rate of enzy- 
matic reactions (because UJ21 represents the rate of ATP 
hydrolysis catalyzed by the CRE motor) [40|. In that 
case a;2i itself would saturate with the increase of ATP 
concentration, instead of increasing linearly with [ATP] . 
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FIG. 10: Normalized MFTT aN~^U plotted against the nor- 
malized motor speed 1^(2/0 for different values of uj2i/a with 
(a) L0i2/a = 0.8, and (b) uji2/a = 0.1, K = 10. 



FIG. 11: Normalized MFTT aN~^td plotted against the nor- 
malized motor speed wfj/o^ for different values of ljju/ol with 
(a) uj2i/a = 0.8, and (b) uj2i/a = 0.1, K = 10. 



IV. CONCLUSION 

In this paper we have studied the process of ATP- 
dependent chromatin remodeling. For simplicity, we have 
considered only a single nucleosome consisting of a ds- 
DNA strand wrapped one and three-fourth turns around 
a cylindrical spool made of histone proteins. We have 
extended Chou's model [23| by assigning two distinct 
"chemical" states to the CRE and postulating a minimal 
mechano-chemical kinetic scheme for capturing the ef- 
fects of ATP hydrolysis explicitly. Our theoretical frame- 
work has been developed exploiting a close analogy with 
the unzipping of a double-stranded DNA by a helicase 
[38|. We have written down the master equations for 
the postulated kinetic scheme. This model of footprint 
traversal by ATP-dependent CRE can be easily inter- 
preted as an implementation of a Brownian ratchet mech- 
anism. From an analytical treatment of this stochastic 
kinetic model, we have derived analytical expression for 
the MFTT of the ATP-dependent CRE. We make explicit 
analytical predictions on the dependence of the MFTT on 
(i) the unhindred speed of the CRE, as well as on (ii) the 
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FIG. 12: Normalized MFTT aN'^td plotted against the ATP 
concentration [ATP] for 0112/0. = 0.1, K = 10.0 and two dif- 
ferent values of L1JI2/0.. 
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concentration of ATP. In principle our theoretical predic- 
tions can be tested by carrying out in-vitro experiments 
with a single nucleosome. 
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